A NOTE ON WIENER-HOPF DETERMINANTS 
AND THE BORODIN-OKOUNKOV IDENTITY 

Estelle L. Basorf] and Yang Chen 

The continuous analogue of a Toeplitz determinant identity for Wiener-Hopf oper- 
ators is proved. An example which arises from random matrix theory is studied and 
an error term for the asymptotics of the determinant is computed. 

1 Wiener-Hopf Determinants 

Recently, a beautiful identity due to Borodin and Okounkov was proved for Toeplitz deter- 
minants which shows how one can write a Toeplitz determinant as a Fredholm determinant. 
In this note we generalize this to the Wiener-Hopf case. The proof in the Wiener-Hopf case 
follows identically with the second one given in |l| . We include it here for completeness sake 
and because the nature of the identity is slightly different in the continuous verses discrete 
convolution setting. 

In the Wiener-Hopf case we begin with a Fredholm determinant on a finite interval and 
then show how this can be written as a Fredholm determinant of an operator defined on L 2 
of a half-line. The point is the second operator has a very "small" kernel and thus higher 
order approximations (as a function of the length of the finite interval) can be found. 

We now state the analogue of the identity and then apply it to a particular case to show 
how error estimates can be computed. In the future we hope to refine the estimates given 
here, apply this identity to other important examples, and also extend it to other operators. 

We consider the Fredholm determinant of the finite Wiener-Hopf operator 



det(/-/C [0 , Q] ), (1.1) 
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where the K, acts on L 2 (0, a) and has kernel K{x — y) with K, given by the Fourier transform 
of a function F, i.e. 

poo 

F(0 = / K{x)e* x dx. (1.2) 



The continuous analogue of the Borodin-Okounkov identity is given, under appropriate 
conditions by the formula 

det(J - K [0M ) = Z e ca det(J - C [a>oo) ), (1.3) 

where 

c:= (1-4) 
C is an operator acting on L 2 [a, oo) with kernel, 

and Z is a certain constant whose value will be defined shortly. Here X is inverse Fourier 
transform of 0(£) , and <f>±(£) are the Wiener- Hopf factors satisfying; 

i-F(0 = <MO<MO, eeR, (1.6) 

and also satisfying the condition that the functions <p± when extended away from R are 
analytic in the upper and lower half-plane respectively. So that this makes sense and our 
proof is valid we require that F is bounded and in L l {R) and that /C is also in £ X (R) and 
satisfies ffj^ \x\ \fC\(x)dx < oo. To achieve an unique factorization, 1 — F must have index 
zero, be bounded away from zero and we assume that <p± are one at ±oo. These conditions 
guarantee that the operator A — I in the proof below is trace class and also that all the 
integrals defined are finite. For details see f|J. 

In this section we include the proof of the identity as promised. In the next section we 
apply it to a particular example that arises in random matrix theory ||. In this example 

vt \ ^sinvr(x-y) 

K,{x - y) := — — -, x, y e R, g > 0. (1.7) 

7rsmhg(x — y) 

We show that as a — > oo, 

det(J - C [a>oo) ) ~ 1 - C(g)e- 2 ^ 1 - e ^\ (1.8) 

where C (g) is a completely determined constant and cos 9 ■= e-* 2 /°, < 9 < tt/2. This is a 
refinement of the classical Szego-Kac-Widom Theorem. Previous attempts at the refinement 
computed the exponential term but not the constant. In principal, as the reader will see, 
higher order terms can also be found using the method outlined in section two. 



Here is a proof of the above identity in the Wiener- Hopf case. As already stated this proof 
follows from []]]. The interested reader should also note another slightly different and very 
elegant proof given in Denote by P a the orthogonal projection of L 2 {0, oo) onto L 2 {0, a), 
Q a = I — P a , and P the orthogonal projection of L 2 (— 00,00) onto L 2 (0, 00). Also define 
W(4>) to be PM^P and W a {(fi) to be P a M^P a where is multiplication by 0. Note that 
by using Fourier transforms it can be proved that W a {<fi) is unitarily equivalent to J — /C[o, a ] 
with <p — 1 — F. It is straight-forward to check that 

P a W(<j>+) = P a W(<f>+) P a , W{<f>.) P a = P a W(</>-) P a , 

and 

w \<t> + )w (cp- 1 ) = J, wtejwtez 1 ) = 1. 

Using the above, we can write Q 

P a W(<j>) P a = P a W(<f> + ) W^Z 1 ) W{cj>) W^Z 1 ) W((f>-) P a 

= P a W{<f) + ) P a W^Z 1 ) W{<j>) W^Z 1 ) Pa W{</>-) P a . 

Now it can be shown that the product of the determinants of P a W((f)±) P a are equal 
to e ca , (see (§], section 10.79). Thus to compute det(J — /C[o, Q ]) we need to consider 

PaW^W^W^Pa- 

Set 

W{<j>^)W{<j>)W{<j>Z}) = A. 

Notice that the determinant of P a AP a equals det {P a AP a + Q a ). Now A is invertible and 
differs from / by a trace class operator ||. Therefore 

det {P a AP a + Qa) = det A det (A _1 P a AP a + A~ l Q a ) 

= det A det (A' 1 (I - Q a )AP a + A _1 Q a ) = det A det (P a - A~ x Q a APa + A~ l Q a ) 

= det A det (P a + A _1 Q Q ) det (/ - Q a AP a ), 

since P a Q a = 0. The determinant of the operator on the right equals one, again since 
PaQa = 0. Moreover 

det {P a + A~ l Q a Q a ) = det (/-(/- A~ l )Q a ) = det (/ - Q a {I - A- l )Q a ). 

We have shown 

det (J - K [QM ) = det A det {I - Q a {I - A- l )Q a ). (1.9) 

2 It is an easy general fact that if tpi € H°° or ijj^ £ then W(i/)tip2) = W (ipi)W '(^2) . In particular 
W(<^>±) are invertible with inverses VF^Ijj 1 ). Recall that H°° consists of all ip £ L°° such that the Fourier 
transform of if) vanishes on the negative real axis. 



It remains to show that this is the same as (1.3). First, A is similar via the invertible operator 
W((j) + ) to W(<f>) W((j)Z l ) W(M)- Therefore 

det A = det W(<j>) W^Z 1 ) W^Z 1 ) = det W(<f>) W(M). (1.10) 

This is a representation of the constant in (1.3). Next 

A' 1 = W ((/>-) W{<j>)~ x W {</>+) = W(<f>-) W^Z 1 ) W{4>Z X ) W((f) + ) (1.11) 

= ^(0_/0+)^(0+/0-)- (1.12) 

Because </>_/</>+ and are reciprocals it follows easily from the algebra properties of 

our operators that the determinant of 

(/ - Q a (I - A-^Qa) 

is the same as the determinant given in the right-hand side of (3) and this completes the 
proof. 

2 An Example 

Next we turn to the example cited in the first section. Consider (1.1) with 

K(x-y):= ^ g>Q (213) 

7rsmng(x — y) 

The Fourier transform is given by, 

F(£) = sinh(7r 2 /(7) 

cosh(7r 2 /g) + cosh(7r£/g) 

It is clear that F satisfies our hypothesis. Writing, 1 — F(£) = exp(ij)(£)), and 

roo rO 

V>(f)=/ Ae m dt+ ipte^dt, (2.15) 

J J -co 

where 

& = [°° ln(l - F{t))e-*p- (2.16) 



yields 

MO 



<M0 



exp(±i$(0), (2.17) 



3 It is known that the constant Z can also be expressed as exp J °° zln(l — F) z ln(l — F)- Z dz, (see 
section 10.79). 



with 

POO 

$(0 = 2% / i) t sin( &)dt. (2.18) 
Jo 

A simple calculation shows that 

C(x,y)= / f(x + z)f(z + y)dz, (2.19) 
Jo 

where 

'MHW^'-^ft-i (2 - 2o) 

To find f(t) we need to compute ip t and $(0- With a : = e -* 2 /\ a calculation gives, 

cos(7rt) - cosh(g6>t/7r) 

V>t = ' . ^ a=:cos#, 0<#<tt/2. 2.21 

tsmh(gt) 

It is a bit difficult to compute $(£) directly, so instead of finding $(£) we attempt to find 
its derivative. Using (2.18) and (2.21) and we find, that 

^ _ . r°° f cos(irt) — cosh(q9t/ir)\ . 

*<«)= 2! i ( ; si nh (gf) j sinKi)tM - (222) 

Next using the above and formula 3.524.5 of J|, gives the integral as a limit as \i — > 1 of the 
sum of C(/i, (1 — /3/j)/2) and £(/i, (1 + /3/j)/2). Then using the definition of the ( function, 
rearranging the sums, and finally letting // — > 1 we obtain 

W0-.g(i + |)-|)--g + «- + ^ 



25 


; 


i (7T — 









Z (7T + 





2^ 




Z (7T — 





2<? 



2 \ 7T y 2gJ \2 

where \I/ is the di-gamma function. Integrating (2.23) with respect to £ with the initial 
condition $(0) = 0, gives, 

exp(i$(f)) 



ri 


(K 1 + J) + 1) 


i n 




r i 


(id +f)-5) 


in 


(Id -!)-|) 



r Q(i + f)-g)r(Ki-f)-g 
r(l (i-f) + |) r(|(i + f) + | 



(2.24) 



We need now to compute (2.20). Note that the integrand vanishes when £ is zero. Further- 
more, by using the asymptotics of the Gamma functions, it is easily seen that the integral 
converges conditionally. Putting s — || and z = exp(—2gx), in (2.20), we find, 

fix) f°° / r(a + s )r(l-a + S ) T{b-s) Y (1 - b - s) \ s ch_ 
2g J-ioo \T(a-s)r(l-a-s)r(l-b + s) T (b + s) J Itxi 1 ' 

where 

a = 




To compute the integral (2.25) we take a sequence of contours consisting of the line segment 
from — %R n to %R n and the semi-circle of radius R n in the left-half plane centered at the 
origin. Here R n = —n + 5, where 5 is any fixed constant satisfying < 5 < 1/4. Now 
computing the residues of the Gamma functions in the left-half plane we obtain after taking 
the limit n — > oo the sum of two hypergeometric functions 4F3; 



m / r (f)r(i-|(f-f))r(i-|(f + f)) 
2s ^ r(i-f)r(i(f + f))r(|(f-f)) 

xe"^ 1 ^)- 4F3 (f^^I;exp(-2^))) 

rH)r(i + |(f-f))r(i + |(f + f)) 
r(i + f)r(|(f -f))r(i(-f -f)) 

xe -K 1+ ^) Q 4 F 3 f ^l 2 '^; exp (-2</a) 

\ 7T ' 5 7T 



(2.26) 



+ 



where 



2 V 71 " £ / 
1 1 /0 m\ 
7i = 1 + o > 72 




As a — > 00, 



where 



tiC [a>oo) ~ C(<?) exp I -2</ I 1 - - j a j , (2.27) 



c , , = , r (f)r(i-|(f-f))r(i-|(f + f) ) 



We have obtained (2.27) from (2.19) and the trace. Note that since < Q/n < 1/2, the 
leading exponential factor, exp(— g(l +0/ir)cx), in the second term of (2.26) will tend to zero 
faster than the first. Furthermore, since 



oo 

4*3 (£::::;£;*) = i + E^ n > \A <i, 



n=0 



where the above series converges absolutely in the unit circle, we see that only the exponential 
factor in the first term of (2.26) will be relevant in the computation of tr£[ a oo ) for large a. 

Here we give an expression for trC^ aoo y where k is any positive integer and with this 
determine the higher order terms for large a. We can compute det(J — £[ aj00 )) by computing 
traces. The justification for this fact follows from our estimate for f(x) when x is large. It 
shows that the function is exponentially small. Hence when x is larger than a it follows that 
the operator £r Qj00 ) has norm smaller than one. 

To compute the traces we note that from the expression of C(x,y) given by (2.19) and 
by shifting the interval of integration, [a, oo), in the trace to [0, oo), 

/*oo /*oo 

trC\ a oo) = ... j f(x! + x 2 + a)...f(x 2 k + X! + a)dx 1 ...dx 2 k- (2.29) 

For fixed x, y > 0, as a — > oo, 

r(l\r(] + l(fL _ «mr(i _ 1(1 . in)) 

/<— ) ~ ^ ( ^ { 44; i( ;': t) ; )1 ^ ». <»•> 

where the second exponentially decaying term in (2.26) has been discarded. Integrating 
(2.29) with (2.30) we find, as a -»• oo, 

tr£f a>oo) ~ C(# exp (-2 5 (l - ^ ko^j . (2.31) 



Therefore, 



oo 4- r /*k 

-lndet(J-£ [a>oo) ) = 

fc=i ft 



oo /^yfe 



E 



-c 



-2 fl (l-|)afc 



Finally, 



det(J - £ [Q;00) ) ~ 1 - C( i7 )e- 2 ^( 1 ^), a ^ oo. (2.32) 



In random matrix theory E(a) := det(7— /C[o,<*]) is the probability that an interval [—a/2, a/2] 
(after a suitable scaling) is free of eigenvalues. Putting g — 0, in (2.1), we have the sme 



kernel, "^g^ . It was shown in || that R(a) := — d '"^^ , satisfies a particular Painleve 
V equation. For a simpler derivation of this and a review of random matrix theory see [[l(J . 
Recently, it was shown in |J using the theory of ly] that i?(a) for g > 0, satisfies a particular 



Painleve VI. In a heuristic perturbative calculation on the Painleve VI both the exponential 
decaying terms were was found |7j], however, the prefactor, C(g), cannot be determined as it 
involves an indefinite integral. It appears that C(g) can only be obtained from determinant 
identity, (1.3). It can also be shown that the computation of (1.1) with the kernel given by 
(2.13) can be reformulated as an equivalent 2x2 matrix Riemann-Hilbert problem, see |j) 
for a description of this technique. Although as a — > oo, such a technique could reproduce 
the classical results of Akhiezer, Hirschman and Szego, namely, the Z factor, it is not at 
all clear whether the higher correction terms given by (2.34) can be found. Final remark: 
the determinant identity, (1.3), does not hold for g — 0. Since in this situation, the Fourier 
transform of the sine kernel is the characteristic function of [— 7r/2, 7r/2] and the Wiener-Hopf 
factorization fails. 
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